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1. Introduction 

The problem of controlling elementary chemical reactions, and especially the advances in 
the field of quantum informatics 1 1 1 stimulate interest towards controlling quantum state of 
microscopic spin systems ||2|. Spin correlated radical pairs are one of the most intriguing 
quantum objects. The quantum state of such a pair can be controlled by applying external 
magnetic fields |4||. Additional application of resonant microwave fields substantially 
widens the possibilities for manipulating the spins. 

Figure[l] shows the scheme of transitions between spin states of a radical ion pair induced 
by hyperfine interaction with magnetic nuclei and resonant microwave field. Usually the set 
of singlet state of the pair \S) with zero total spin and three triplet states |T+), \Tq) and |T_) 
with total spin 1 and its projection on the z axis equal to +1, and —1, respectively, is taken 
is the basis. In strong applied magnetic field the degeneracy of the triplet states is lifted by 
Zeeman interaction. In the absence of microwave field and relaxation the populations of |T+) 
and |T_) states, as well as the sum of the populations of |To) and \S) states remain constant. 
In the presence of hyperfine interaction the states \Tq) and \S) are not stationary, and their 
populations periodically change with time. This phenomenon is referred to as quantum beats 
and is observed in experiment ||5j|6|. Quantum beats can also be induced by differences in the 
(yf-values of the pair partners fT). 

A resonant microwave field induces transitions between the triplet states, which leads to 
changes in the population of singlet state of the pair and can be used to optically detect ESR 
spectra of radical ions in liquids |l8l|9|. If the rate of transitions between 15) and |To) states 
is higher than between the triplet states, periodic changes of the singlet state population with 
frequency proportional to amplitude of the microwave field are observed, which is referred to 
as quantum oscillations ifTOlfTTllTll . In the opposite limiting case of strong microwave field, 
when the rate of transitions between the triplet states exceeds the rate of transitions between 
\S) and I To), singlet-triplet transitions are slowed down, and the so-called spin locking is 
observed. As was demonstrated [13| for pairs formed in singlet initial spin state a strong 
microwave field only slows down singlet-triplet transitions but does not completely block 
them, and the population of \Tq) state always remains low. This was suggested as a possible 
route to control spin state of such a pair. If the microwave field is switched off at the moment 
when the population of the singlet state is minimal, most pairs will be trapped in the |T+) and 
|T_) states. The only way to get into the singlet state for them is spin-lattice relaxation. A 
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Figure 1. Spin dynamics of a radical ion pair 
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similar idea was exploited to experimentally substantially prolong lifetime of the radical pair 
in photosynthetic reaction center [141 by transferring it into its |T+) and |T_) states. Another 
route to control spin state of the pair was found theoretically in the work lITSll . If a strong 
resonant microwave field is rapidly switched on at the moment when spin system is in its 
|To) state, the pair is completely locked in triplet state and will never become singlet through 
dynamic transitions. A possibility of controlling spin state of the pair with very short pulses 
of very strong microwave field has also been treated theoretically 1 16|. The theory shows that 
quantum beats in radical pair can be controlled with two (for singlet-born pairs) or just one 
(for triplet-bom pairs) short microwave pulses. 

The already cited works lfT3lfT5ll treated spin dynamics of a radical pair in static magnetic 
and resonant microwave fields neglecting relaxation. A method for solving the problem of 
spin dynamics taking into account relaxation was developed in work [12J, and approximate 
solutions in the two limiting cases of large and small splitting were obtained. A method for 
numerical solution of this problem taking into account relaxation and ion-molecular charge 
transfer reaction was developed in work ifTTll . In the present contribution we provide a general 
analytical solution for the problem of spin evolution of a radical ion pair taking into account 
spin relaxation and use it to analyse the characteristic features of time evolution of the pair 
spin state in microwave field. 

2. Analytical solution 

2.1. Model 

Let us consider what happens when a non-polar solution is irradiated by ionizing radiation. 
As an example we shall use liquid hydrocarbons (RH), in which ionizing irradiation in the 
presence of electron (A) and hole acceptors (D) initiates the processes which in turn lead 
to formation of several types of singlet spin correlated pairs (RH-+/e~), (RH-"^/A-~), 
(D-+/e~) and (D-+/A-^) and their subsequent recombination [fT8lfT9l. Singlet excited 
molecules (^A* and ^D*) produced by recombination emit a quantum of light that is detected 
by experimental setup. The intensity of the detected luminescence is thus proportional to 
probability of forming singlet excited molecules upon recombination. The effects of magnetic 
fields on recombination luminescence have as their origin spin evolution in geminate pairs and 
dependence of the yield of luminescence on the multiplicity of excited molecule formed by 
recombination. 

An important advantage of non-polar solutions is that the initial inter-partner distance in 
the pairs after ionization (normally ~ 5 6 nm) is substantially lower than Onsager radius 
(~30 nm for alkanes at room temperature). The overwhelming majority of the pairs thus 
recombine as geminate pairs, i.e., with their sibling counter ion. Upon ionization of molecule 
the spin of the ejected electron does not change its state. Since in the molecule the spins were 
paired, the initial spin state of the geminate pair is always singlet. Charge transfer to acceptors 
and charge recombination also do not change spin states. This means that the multiplicity 
of the excited molecule formed by recombination is determined by the multiplicity of the 



Time evolution of spin state of radical ion pair in microwave field 



4 



recombining pair immediately before recombination. The multiplicity of the excited molecule 
is thus determined by spin evolution of the geminate radical ion pair between the moments 
of ionization and recombination. Static or oscillating magnetic fields can substantially affect 
this evolution. And finally, as excited molecules normally emit light from their singlet excited 
state, any effect of magnetic fields on the spin state of the pair is directly reflected in the 
intensity of recombination luminescence. 

In liquid non-polar solution the partners of the pair formed by irradiation spend almost 
the entire period of time from ionization and up to recombination at substantial distances 
from each other (tens of nanometers), as their approach to distances about 1-2 nm leads to 
practically instant recombination. Thus dipole-dipole and exchange interactions between the 
pair partners can be safely neglected. Spin evolution of the pair in this case is driven by 
interaction of electron spin with spins of nuclei (hyperfine interaction, HFI), with external 
magnetic fields, and by spin relaxation processes. 

Spin Hamiltonian for radical pair consisting of radicals A and D in high static magnetic 
field Bq and microwave filed with magnetic component Bi{t) can be written as follows: 

H = g/3[Bo + B,mS^ + S^) + Sf ^ afl^ + Sf ^ af Jg, (1) 

i j 

where and are spin operators for electrons, and If- and /j^ are spin operators for 
magnetic nuclei in radicals A and D respectively. The sum runs over the nuclei of radicals. 
The first term describes Zeeman interaction of electron spins with external magnetic field, the 
second and third terms correspond to isotropic hyperfine interaction of electron spins with 
nuclear spins in the two radicals. The z axis is aligned along the direction of the field Bq. For 
simplicity we assume that the partners of the pair have identical g-values (gA = go = g)- 

As was shown in work 11121 in the presence of spin relaxation time evolution of the spin 
operators for the partners of the radical pair S^{t) and S^it) is describe with a system of 
equations in the frame rotating at the frequency of the applied microwave field that are similar 
to Bloch equations 1.201 l2T1l . For operator S"^{t) the system is written as: 

dSf/dt = - Auj^Sf - Sf/Tf, (2a) 
dS^/dt = Ato^Sf - uiSf - S^/Tf, {2b) 
dSf/dt = uo^S^ - Sf/Tf, (2c) 

where uji = gf3Bi/h, Tf and Tf are the times of spin-lattice and spin-spin relaxation. 
The z axis is aligned in the direction of the static magnetic field, and the x axis points 
along the Bi field. Au;"^ = cu"^ — luq is the detuning of the hyperfine component (HFC) 
of radical A from resonance. Here luq is the frequency of the applied microwave field, and 
uj^ = {gPBo + '^ciflf)/h is the resonance frequency for radical A. Time evolution for 
partner S^(t) is described by similar equations after substituting D for A. 

In the absence of interactions between the spins of the pair spin state of the pair can be 
described using operators of projection ^/.(t) on an arbitrary state \ip) written through 5^(t) 
and S^{ty. 

w4t)=Tr[P4t)m], (3) 
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2.2. Solution of the system of operator equations 



Equations (I2al) -(l2cl) are a system of homogeneous linear differential equations, and its 
solution can be represented as an expansion over a complete orthonormal basis. In our case 
such a basis is conveniently given by Pauli matrices o"^, a^^, and a^. 

St'''{t) = Y,^''it)Sm, hk = x,y,z, (4) 

k 

where 5^(0) = |(Ti, i = x,y, z. 

The system (Oal) -(Ocl) was solved by Laplace transform and its inversion along the lines 
of the work |22|. A substantial difference from \t22i is that in our case the equations are solved 
for operators, but this does not change the basic approach. 

Laplace transform of the original equations yields the following system: 

-AujS., + ip+ l/T2)Sy + UJlS, = Sy{0), (5b) 

-uJi§y + {p+l/T,)§, =4(0), (5c) 

where the tilde sign denotes Laplace transforms of the corresponding spin operators, p is the 
Laplace parameter. This system is a system of linear algebraic equations with three unknowns, 
which has a unique solution when its determinant is nonzero 

A{p) = {p+ im'ip + l/Ti) + Auj\p + 1/Ti) + Lulip + I/T2). (6) 

The solution of the system (I5at-(l5cl) has the following form: 

/ iP+T,)iP+jr)+^! -ip+jr)Au AUU, \ 

{p+^)Au -{p+^)u, 5(0) (7) 

V ^^^1 (p + i)^i + {p + ) 

To take inverse Laplace transform from expressions dTJ) we need to factorize the 
determinant A(p). The determinant is a cubic polynomial with respect to p, and its roots 
are given by Cardano formula. Depending on whether the discriminant of the equation 

e = gV4 + sV27 (8) 

is positive, negative or zero (two cases of degeneration), four types of solution are possible. 
We introduced the following notation: 

q = 26V27 - bc/3 + d, s = -^^/S + c, b = 2/T2 + 1/Ti, (9) 
c = l/r| + 2/(TiT2) + Au^ + ul, d = l/(T|Ti) + Au^Ti + U0I/T2. (10) 
As can be seen from equations dTJ), in all cases the following relations hold true: 

A^yx{t) = -Ky{t), X,^{t) = Xx,{t), X,y{t) = -Ky,{t). (11) 



1 



A{p) 
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In the case ^ > the polynomial A(p) has one real and two complex conjugate roots 
(the roots of A(p) are —pi and —p2 ± ips): A(p) = {p + Pi)[ip + P2Y + pI], where 
p, = b/3-{y++y^), P2 = 6/3+(y++Z/-)/2, ^3 = VS{y+-y-)/2, y± = {-q/2 ± ^f^' . 
This case corresponds to an oscillatory solution: 

A,fc(t) = Ake-P'' + B^ke-P^' cosip^t) + -^^6-^^* sin(p3t), (12) 

P-i 



The coefficients were found by procedure similar to the method used in the work II22I . The 
expressions for the coefficients are given in [Appendix A 



The solutions for species type A and type D are similar and differ only in the 
corresponding indices (A or D) for parameters Acu, Ti and T2. 

[Appendix Aj does not provide expressions for coefficients Ayx{t), Azx{t) and Azy{t), as 
they can be easily obtained from formulae dTTT) . 

In the case of negative discriminant (^ < 0) the polynomial A(p) has three real roots: 
A{p) = {p + Pi){p + P2){p + Ps), where Pi = 6/3 + 2pcos(0/3), p2 = b/3- 2pcos(7r/3 - 
(/,/3), p3 = 6/3 - 2pcos(7r/3 + (j)/3), p = ±^-s/3, cos(0) = g/(2p3) and the sign of p 
must coincide with the sign of q. Similar to the first case we obtain: 

A,fc(t) = Cl,e-''^\ps - P2) + Cle-^-\p, - p,) + Cle~^^^\p2 - Pi). (13) 

This case corresponds to a overdamped solution (all pi > 0). 

The expressions for the coefficients C'^ are given in [Appendix A Coefficients Ay^^t), 
Azx{t) and Azy{t) can be found using formulae (fTTTl . 

For ^ = the determinant A(p) has three real roots, two of which coincide: A(p) = 
(p + Pl)(p + P2)^ where pi = 6/3 - 2(-g/2)^/^, p2 = 6/3 + (-g/2)^/3. In this case: 

A.fc(t) = Afce-P^* + Bike-P^' + Cfcte"^^*. (14) 

For ^ = and q = the determinant A{p) has three coinciding real roots: 

A{p) = (p + pi)3, wherepi = 6/3 = (2/T2 + l/Ti)/3. In this case: 

A,fc(t) = AkB-P^' + B^kte^P^' + Cj^e-^^'. (15) 

This case (double degeneration) occurs when lui = 2\plAuj and (I/T2 — 1/Ti) = 3\f3Aijj. 

The degenerate cases ([T4b and dTSt never occur in actual calculations, so we do not 
provide expressions for their coefficients. 

3. Some peculiarities of spin dymanics 



A convenient basis for spin states of the radical pair is given by Bell functions 

1 

1 

1 
71 



\s) = \r) = - \m)\ (16a) 

T.) = |0-) = ^(IT)IT)-|i)|i)), (16Z7) 
T,) = |0+) = 4(lt)IT) + U)U)), (16c) 



Time evolution of spin state of radical ion pair in microwave field 

1 



\T.) = \r) 



V2 



(IT)U) + U)IT)), 



(16 J) 



where | |) and | J,) are eigenstates of radical with spin projection on z axis equal to +1/2 and 
—1/2 respectively. In the product the left function always corresponds to radical A and the 
right function to radical D. The singlet function \S) has zero total spin, and the three triplet 
functions |Tj.), \Ty) and IT^) have total spin equal to 1 and zero projection on axes x, y and z 
respectively. 

Specialists in spin chemistry are more used to working in the basis of the singlet function 
\S) and the three triplet functions \Tq), |T+) and |T_) with defined projection of the total spin 
on the z axis equal to 0, +1 and —1 respectively (it is this notation that was used in Figure[l]). 
As equations (I16Z7|) - (I16JI) show, the Bell functions are expressed from them as follows: 

1 , ,„ 1 



|T,) = |To), |T,) 



V2 



(|T+)-|T_)), |T,) 



71' 



\T+) + \T^)). 



The expressions for spin projection operators Pa in the basis \S), \Tx), \Ty), \Tz) is given 
in [Appendix C 



m 



(17) 



Let us introduce notation for density matrix written in this basis at time zero: 

^ Pss Psx Psy Psz ^ 
Pxs Pxx Pxy Pxz 
Pys Pyx Pyy Pyz 
y Pzs Pzx Pzy Pzz 

In our calculations we shall also use a particular case when initial state of the pair at zero 
time is given by a wave function that is a linear combination of singlet \S) and triplet \Tz) 
states: 



|^)(0) = cos^|5) + sin^e^'^|T,) 

In this case pss = cos^ 9, p^z = sin^ 9, psz 
equal to zero. 



I sin 2^6" 



(18) 

i sin26'e*'^, all other pij are 



3.1. Spin dynamics in the absence of spin relaxation 

In the absence of spin relaxation, i.e. for 1/Ti = I/T2 = 0, the calculations become 
substantially simpler, which allows to get a number of important results analytically. In this 
an oscillatory solution always occurs : 

Aifc(t) = Aik + Bit cos(u;t) + —Cik sin(u;t), (19) 

where uj = a/cUi + Auj"^. 

3.1.1. Spin dynamics in the absence of microwave pumping (uJi = 0) As it follows from 
expressions (IC.1I) - (IC.4I) . 

Pss + Pzz = \S){S\ + |T,)(T,| = 1 - 2S^it)S^{t), (20) 
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P.. + P» = |r.>{r.| + |r,)(r,| = i + 2S^WSfW. (21) 

For each radical we can write down: 

SM = (1 - cosM)) SM + — sinM)4(0) + + 4 cosM)) 5,(0). (22) 

For cji = and considering ( IC.5I )-( IC.13I) we have: 

Tr[5,^(t)5f (t)p(O)] = i(-p,, + p., + p,, - p,,), (23) 

where we used the notation of (flTl) . 

From (l20b . (l23t and from the normalization condition 

Pss + Pxx + Pyy + Pzz = 1 (24) 

it follows that for ui = the sum of populations of the singlet state pssit) and teh state IT,) 
Pzzit) does not depend on time: 

Pss(t) + Pzz{t) = Pss + Pzz- (25) 

This is a quite natural result, since in the absence of microwave pumping and spin relaxation 
the only possible transitions in the system are between the states 15) and IT,) due to hyperfine 
interaction or differences in the values of the pair partners. The two other triplet states 
|T+) and |T_) (that can be linearly combined to produce the \Tx) and \Ty) states) in these 
conditions are stationary, and their population does not change in time due to substantial 
Zeeman splitting. 

3.1.2. Spin dynamics in a strong microwave field (Au/uji 0) As it follows from 
expressions (IC.1I) - (IC.4I) . 

Pss + p.. = \S){S\ + \T,){T^\ = i - 2S^it)S^it), (26) 

Pyy + Pzz = \Ty){Ty\ + |T,)(T,| = i + 2S^{t)S^{t). (27) 
Since for each radical 

SAt) = f4 + ^ cosicot)] 5,(0) - — sin(a;t)4(0) + (1 - cos(cut)) 5,(0), (28) 

for /\uj /uji — > considering (IC.5I) - (IC.13I) we have: 

rr[5f (t)5f (t)p(O)] ^(-p,, - p,, + Pyy + Pzz). (29) 

From (l26b . (l29b and (|24h we obtain that in the limit of strong microwave field a new 
invariant appears. The sum of populations of the singlet |5) and the triplet IT,) states does 
not depend on time: 

Pss{t) + Pxx{t) ~ Pss + Pxx- (30) 

The equality becomes the more accurate the smaller is the ratio Acj/cui. 
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3.1.3. Spin dynamics in microwave fields of finite strength The situation in arbitrary 
microwave fields is more complex than in the limiting cases discussed above, so we shall 
consider the general properties of the obtained solutions using several specific examples as an 
aid. 

Figure Eta) shows time dependencies of the sum of populations of the singlet \S) and 
the triplet \Tz) states of the system at different magnitudes of Bi. The line of one radical is 
precisely in resonance with the applied microwave field, and the line of the other partner is 
detuned by 0.3 mT. Initial state of the pair was taken to be singlet, i.e. pss = 1, all other 
Pij = 0. For Figure 121 and all next figures are true ABa,d = hAcu^'^ /{gP). As the figure 
shows, in this case even a weak microwave field is sufficient to completely destroy invariant 
(l25t . The reason for this is transitions between jT^) and \Ty) states induced by resonant 
microwave field. As a consequence the value of pssit) + Pzz{t) periodically changes from 
zero to one. The population of the singlet state also changes periodically (see Figure |71). This 
phenomenon is referred to as quantum oscillations and has been observed experimentally 
lfT0lfTnfT2t and studied theoretically lfT3lfT2t . The changes in the frequency of oscillations 
with increasing amplitude of the microwave field are not monotonous. The frequency first 
increases but then, as Bi grows substantially larger than the splitting between the lines of the 
two radicals, falls down. The decrease of the frequency is related to the effect of spin locking. 

Figure Etb) shows time dependences Pss{t) + Pzz{t) in the situation of strong detuning 
of the microwave frequency from resonance. Time scale is expanded three five times with 
respect to Figure IJa). As is clear from the figure, in this case for not too large amplitudes of 
the microwave field the frequency of oscillations substantially increases as compared to the 
case of resonance with accompanying decrease in the amplitude of the oscillations. 

This effect can be observed experimentally. Figure |3] shows experimental results on 
monitoring kinetics of recombination fluorescence from liquid alkane solutions in the absence 
and in the presence of microwave field. The experimental technique that was used to get 
the presented traces was described in details in work ifTlll . The processes initiated in the 
sample under the action of a short pulse of X-irradiation are given by scheme (??). Dodecane 




Time (^is) Time (us) 

Figure 2. Time dependence of the sum of populations of the singlet \S) and the triplet \Tz) 
states. ABa=0, ASd=0.3 mT (a); AB^=1 mT, AB£,= l-3 mT (b). The value of Bi is 
indicated next to its corresponding curve 
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was used as the solvent (RH), hexafluorobensene CqFq served as the electron acceptor 
(A), and deuterated jpara-terphenyl PTP — di^ as the hole acceptor (D). ESR spectrum of 
PTP — c/i4 radical cation is a single inhomogeneously broadened line with width about 0. 1 mT. 
Radical anion of CqFq has a wide resolved spectrum with hyperfine coupling constant of 
about 13.5 mT with six equivalent fluorine nuclei f23\. Time dependence of the intensity 
of recombination fluorescence I(t) was sampled using single photon counting technique. 
The observed light in this case mostly comes from excited PTP — di^ molecules. Since 
recombination kinetics f(t) does not depend on the multiplicity of the pair, and fluorescence 
intensity is proportional to singlet state population of the pair at the moment of recombination 
Pssit), we can write that I{t) ^ f{t)pss{t). The equality becomes more accurate as 
fluorescence time of the luminophore shortens (about 1 ns for PTP — du) and time resolution 
of the experimental setup improves (about 3 ns in our case). In these conditions time resolved 
microwave field effect is given by the following expression: 

ImUt)/Ioff{t) - 1 ^ pTs'"it)/p"/J{t) - 1, (31) 

where Imwif) and are luminescence kinetics and time dependence of the singlet state 

population in the presence of microwave field, and Ioff{t) and p°sl^{t) are the same functions 
in the absence of the field. 

Figure |3] shows three curves. Curve 1 corresponds to conditions of microwave field 
resonance for radical cation of PTP — c?i4, curves 2 and 3 — to detuning from resonance by 
0.3 mT and 1.3 mT respectively. As can be seen from the figure, the experimental curves 
are modulated with oscillations. The frequency of oscillations increases with increasing 
the detuning, and their amplitude falls down in qualitative agreement with the results of 
theoretical calculations shown in Figure |2l Because of rather large HFI constants in radical 
anion of CgFg the populations of \S) and [F^) states can be considered equal at any time, and 
Pssif) = ^[pssit] + Pzzit)]- Exprimental curves are noisy as they were not passed through any 




50 100 150 200 

Time (ns) 



Figure 3. Experimentally observed microwave field effect in the kinetics of recombination 
fluorescence from dodecane solution of 10^* M PTP — du and 10~^ M CgFg in the 
conditions of resonance (1) and with detuning of 0.3 mT (2) and 1.3 mT (3). i?i = 1 mT. 
For convenience the curves 2 and 3 are shifted upwards by 0.05 and 0. 1 respectively 



Time evolution of spin state of radical ion pair in microwave field 



11 




0,0 0,1 0,2 0,3 0,4 0,5 

Time (^s) 



0,0 0,1 0,2 0,3 0,4 0,5 

Time (us) 



Figure 4. Time dependence of the sum of populations of the singlet \S) and triplet |T^) 
states for different amplitudes of microwave field. The magnitude of Bi is indicated next to 
its corresponding curve. ABa=0, Ai?/)=0.3 mT (a); ABa=1 mT, Ai?£)=1.3 mT (b) 



smoothing routines, and the figure shows that noise level increases with time. As was shown 
in work lfT2l . it is proportional to t^^^. 

As Figure Etb) shows for large magnitudes of Bi in the case of detuning from 
resonance the frequency of oscillations also decreases, but this happens at larger amplitudes 
of microwave field than in the conditions of resonance and is accompanied with increasing 
amplitude of the oscillations. 

Figure lit a) shows time dependences of the sum of populations of the singlet \S) and the 
triplet \Tx) states at several amplitudes of microwave field. In the calculations it was assumed 
that one radical is precisely in resonance with the applied microwave field, and the line of the 
second radical is detuned by 0.3 mT. All calculations were performed for singlet initial state. 
Figure Efb) shows the results of calculations similar to those shown in Figure Ufa), but with 
detuning from resonance by 1 mT. 

As FigurelU shows, the behaviour of these curves qualitatively differs from the curves of 
Figure |2l An increase in the amplitude of microwave field leads to increased frequency and 
decreased amplitude of the oscillations. Such a behaviour of the calculated curves reflects the 
existence of the high field invariant (l30b . It can be noted that the accuracy of conserving the 
value of Pss{t) + pxx{t) at small splitting is rather high already for microwave field amplitudes 
easily accessible in experiment. The magnitudes of Bi amplitudes required for conserving the 
high field invariant increase in the case of detuning from resonance (or increasing the value 
of Alu). 

For singlet initial state of the pair and sufficiently large magnitude of Bi the sum of 
populations pss{t) + Pxx{t) oscillate close to one, as can be seen from FigurelU which is in 
good agreement with equation (l30b . In a more general case dTSt . when initial singlet state 
population Pss(O) is equal to cos^ 6, the sum of populations Pss{t) + Pxx{t) should oscillate 
close to this value cos^ 6. The correctness of this conclusion is demonstrated by the results 
of calculations given in Figure |5l which shows time dependences of the sum of populations 
of the singlet \S) and the triplet \Tx) states of the system for different initial conditions (fTSt . 
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Figure 5. Time dependence of the sum of populations of the singlet \S) and the triplet \Tx) 
states for different initial conditions given by equation (I18> . The values of 9 are indicated next 
to their corresponding curves. </> = tt/2 (a) and (b). ABa=0, A_Bd=0.3 mT. Bi=2 mT 



The curves that are shown in this figure correspond to several sets of values of 6 and 0. As 
can be clearly seen, the curves in fact oscillate close to the level of cos^ 9. The frequency of 
these oscillations is indeed close to cui, and their amplitude depends on (p. 

3.2. On possibility to control spin state of the radical ion pair 

The existence of invariants (l25t and (l30b significantly simplifies analysis of spin evolution of 
radical ion pairs in the absence of microwave field and in strong microwave fields. This opens 
important possibilities for controlled impact on the pair to manipulate its spin state. 

As was demonstrated in [|15J . the existence of the high field invariant (l30b allows in 
principle to perform "true spin locking". Upon rapid application of a strong microwave 
field at the moment when the pair is in its triplet state the system is "locked" in triplet 
state and cannot get into singlet state. At the same time for singlet initial state the system 
cannot be completely "locked", and singlet state population oscillates at a frequency close 
to Auj\^/2uji, where Aload = — where and are resonance frequencies for 
radicals A and D respectively lfT3l . This effect is made absolutely clear by the existence of 
the high field invariant (l30b . For singlet initial state the sum of populations Pss{t) + Pxx{t) 
is equal to 1 all the time, and the populations pss{t) and pxx{t) change in antiphase with 
frequency /\uj\j^/2uji. The population of the state \Tz) always remains very close to zero 
11131 [T5I . From the invariant (BUI) it follows that the population of the state \Ty) also remains 
zero. For triplet initial state \Tz) the sum of populations pss{t) and pxx{t) remains zero all the 
time, which means that populations Pss{t) and Pxx{t) remain zero. 

In radiation processes initial state of the pair is always singlet. However the pair can 
be transferred into triplet state jT^) due to existence of invariant (l25t . In the absence of 
microwave field only the states \S) and \Tz) get populated. If the frequencies of transitions 
between them Auad for different pairs in the sample are multiples of each other, as it happens 
when there is an odd number of equivalent magnetic nuclei in the radical ion, there are 
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I 1 I 1 I 1 I I I 1 I 1 

0,0 0,1 0,2 0,3 0,4 0,5 0,6 

Time (us) 



Figure 6. Time dependence of the sum of populations of the singlet \S) and the triplet \Tx) 
states. Initial state is \T^). Bi=l mT, T{^=Tf =10 /is, ABa=0, ABd=03 mT, T^=Tf = T2. 
Numerical values of T2 are indicated next to the corresponding curves 




Time (^.s) Time (^s) 

Figure 7. Time dependence of the populations of the singlet state \S) (1), triplet state \Tz) 
(2) and the sum of populations of the triplet states {Tx) and \Ty) (3) for = T[' = I s, 
= Tf ^1 /IS. ABa=0, ABd = 0.3 mT. Bi=10 ^iT (a) and 1 mT (b) 

moments when all pairs are in their \T^) states. If a strong microwave field is applied at 
this time, the effect of "true spin locking" can be observed. 

3.3. The effect of relaxation 

The above considerations are valid only in the absence of spin relaxation, which destroys spin 
invariants. The zero microwave field invariant (l25t is apparently destroyed by spin-lattice 
relaxation. Phase relaxation does not affect this invariant. This means that expression (l25t 
holds true only for times substantially smaller than Ti. The effect of relaxation on the high 
field invariant (l30b is much less obvious. Figure|6l shows transformations of time dependence 
for the sum of populations of the singlet | S) and the triplet \Tx) states of the pair with changing 
phase relaxation time T2. In the absence of relaxation because of invariant (l30t the total 
population of pss{t) and pxx{t) should oscillate close to zero. However the figure shows that 
this is no longer the case. The total population increases with time, and the time scale of 
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this increase is close to the phase relaxation time T2. We conclude it is phase relaxation that 
causes the destruction of the high field invariant (l30b . 

Figure IT] shows time dependences of the populations of different states in the cases of 
weak and strong microwave field in the presence of phase relaxation. 

The time of spin-lattice relaxation Ti in this model calculation is taken to be very large, 
so it has practically no effect on the obtained curves. 

As Figure |7] shows, time profile of the populations of different states are qualitatively 
different upon application of weak and strong microwave fields. In a weak field pss{t) 
and Pzz{t) oscillate in antiphase with frequency close to /S.ujad- The amplitude of these 
oscillations drops because of phase relaxation, and pss{t) and pzz{t) asymptotically approach 
each other. The sum of populations Pxx{t) + Pyy{^) in this case is a more smooth curve, which 
is periodically modulated with frequency close to uji. 

In a strong microwave field we see the manifestations of spin locking. Singlet state 
population Pss(t) oscillates at a frequency about /\uj similar to the weak field case, but the 
amplitude of these oscillations is rather small. Oscillations of substantial amplitude occur 
at a frequency about ^uj\j^/2uji. This is also the frequency, at which changes the sum of 
populations pxx{t) + Pyy{t)- The population of the triplet state \Tz) increases with the time of 
phase relaxation T2. 

4. Conclusions 

First of all we note that all these results were obtained in the case of fixed nuclear 
configurations in both partners of the pair. In the presence of the reaction of ion-molecular 
charge transfer this condition is no longer observed, but this apparently does not change the 
zero microwave field invariant (l25t . The high microwave field invariant (l30t will also remain 
valid if the condition Auj/uji ^ 1 is not violated in the course of the reaction. 

Let us formulate the main rules for spin evolution of a radical pair that follow from the 
existence of invariants (l25t and (l30b . 

(i) In the absence of microwave pumping the sum of populations of the states \S) and \Tz) 
(or \S) and |To) in different notation) remains constant at times shorter than Ti. The 
populations of these states oscillate at a frequency c<j^ — cu^, where c<j^ and cu^ are 
resonance frequencies for radicals A and D respectively. It also naturally follows that 
the sum of populations of the states \Tx) and \Ty) (or |T+) and |T_)) does not change 
as well. Invariant (l25t is destroyed by spin-lattice relaxation. Even weak resonant 
microwave field also destroys invariant (l25l) . The sum of populations pssit) + Pzz(t) 
starts oscillating at a frequency close to tui (quantum oscillations induced by microwave 
field). 

(ii) In a strong resonant microwave field remain constant the sum of populations of the states 
\S) and \Tx), as well as \Ty) and \Tz). The consequence of this is the effect of spin 
locking, the high field invariant (l30b is rather accurately conserved already for not very 
large magnitudes of Bi quite accessible in experiment. The magnitude of Bi required 
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to maintain the invariant (BUI) increases with detuning from resonance. The high field 
invariant (l30b is destroyed by phase relaxation. 

The two invariants that we found suggest a way of controlling transitions between spin 
states of a radical pair. One possible variant of controlling quantum spin state would be a 
rapid switching of a strong microwave field on and off. In this case the pair will alternatively 
be in the conditions when one of the invariants is conserved. Another possible control option 
may be a rapid change in the phase of the microwave field by 7r/2. In this case the axes x and 
y are exchanged, which correspondingly changes the populations of the IT^.) and \Ty) states. 
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Appendix A. Coefficients in the expressions for different types of solutions of equations 
for spin dynamics 

Case i> Q (for oscillatory solution) 

= (a/5 + ujDh, Ayy = al3h, A,, = {p^ + Acj^)/^, 
Axy = -aAu/j, A^^ = uJiAuj/-f, Ay^ = -uji(3/-f, 

R — 1 — 4 R — ^ — A U — 1_4- 

^xx ^zx) ^yy yy^ ^zz ^zz) 

R— — 4 R— _4 u — _ A . 

^xy ^xj/; ^xz ■'^xzt ^yz yz) 

Cxx = {al3 + ul){pi -V2)h + (a + + -P2), 

Cyy = aj3{pi -P2)h+ {a + 13 + Pi -P2), 

Cz = + Auj^){pi - P2)h + (2/? + Pl- P2), 

Cjcy = -ALj[a{pi - P2) + 7]/7, C^i.^ = ujiAlj{pi - P2)h, 

Cy, = -^1(7 + -p2))/7, 

where 

a = 1/Ti -pi, (3= I/T2 - pi, 7 = (P2 - Pif + pi- 
Case ^ < (for overdamped solution) 

CL = -{Pi + ^)Pi + ^1 + 7f^}, 

V Ti T2 T1T2 



.1 ^1 1^ 1 



Cyy = -APi - (— + —)Pi + 7F7^}, 
1.0 2 .0 1 



2 



ciy = —(p. - cL = -Auu,, cl = ^(p, - 1; 

^ Tj Ti T] " T] T2 

where / = 1, 2, 3, i] = (pi - P2){p2 - P3)iP3 - Pi)- 
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Appendix B. The values of parameters that lead to certain types of solution for 
equations of spin dynamics 

To analyse the different cases of spin dynamics it would be very useful to have simple rules 
for determining the type of solution for the equations with the given set of parameters. As it 
turned out, these rules can be easily formulated, which will be done in this section. 

As has been shown earlier, the exact analytical solution for the equations of spin 
dynamics can belong to one of four types: oscillatory, overdamped, and two degenerate, 
depending in the first place on the value of ^. 

To simplify further treatment let us introduce the new parameter: 

z = I/T2 - 1/Ti. 

Since T2 ^ Ti in all cases, z ^ 0. 

It turns out that relaxation times Ti and T2 enter the expressions for q, s and ^ (l8t- (fT0b 
only as the z parameter: 

1 9 1 

? = 3K - )^ - > * = ~3^ +^1' 

^ = -i- [AAuj^z'^ + (8Acj^ - 20Alj^ujI - uo^)z^ + A{Auo'^ + uolf] . (B.l) 

This means that the type of the sought solution depends only on the differences in the rates of 
spin-spin and spin-lattice relaxation, but not on the absolute values of the two rates. 

As can be seen from expression (IB. II) . ^ is a quadratic function with respect to 2^, with 
branches of the parabola going upwards. This means that in the absence of real roots ^ is 
positive for all values oi z. If the roots exist ^ is negative for values of z'^ in-between the roots 
and is positive for z'^ outside this range. 

From ,^ = it follows: 



"1,2 



(ut + 2QAuj^ujI - 8Acj^ ± A LfK2 - 8Acj2)3 ) . (b.2) 
q/Xuj^ \ J 

The solutions exist for uji ^ 2\/2/\uj. It can be easily shown that both roots are positive. 

From this we conclude that ^ > either for uoi < 2^/2Auj, or for ui > 2\f2Auj when 
z < Z2 ox z > zi {zx corresponds to choice of in (IB. 21) . and Z2 to the choice of "— " 
respectively). In this case we obtain a solution of the oscillatory type. 

^ < for Co"! > 2\f2Aijj in the range Z2 < z < z\. This set of parameters produces a 
solution of the overdamped type. 

For ^ = 0, when either z = or 2; = ^2 we have a solution of the degenerate type. 

Finally, when simultaneously ^ = and g = the case of double degeneracy takes place. 
It can be easily shown that this occurs for cji = 2^/2 Alu and z = Sy/SAu. 

Especially interesting is the behaviour of the solution in the vicinity of resonance for 
one of the radical ions, that is, when Auj /uji 0. Small parameter expansion of (IB. 21) then 
yields: 

Zi ^ ujI/{2Auj) ^00, Z2 ~ Scui. 

Thus close to the resonance line an oscillatory solution will be obtained for coi > \{1/T2 — 
1/Ti), and an overdamped solution for c<Ji < |(1/T'2 — l/Ti). 
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Below we express projection operators of the basis states in the representation \ S), |Tx), \Ty), 
\Tz) in terms of spin operators of the pair partners: 

1 

4 

_ 1 
~ 4 



\S){S\- 



oA oD oA 

X X y y 



qA qD 



qA qD I qA qD I Qji Qi. 



qD 



yy 



\rp \lrp I qAoU qAqU I qAqD 

~ ^ "^y "^y '^'^z'^z 



(C.l) 
(C.2) 
(C.3) 



Pzz 



1 
4 



Irp \/rp I ^1 qA QJJ _|_ QAqJJ qA QJJ 

l-Lz/X-Lz] — ~ ~r Jx 'Jx "T "-"y ^S/ ^2 • 



(C.4) 



Let us express the trace of the product of spin operators by initial density matrix in the 
notation of (fTTl) : 

1 , 

^ Pss Pxx ~l~ Pyy ~l~ Pzzji 



Tr[S^{0)S!'{0)p{0 



y y^J-^y 



Tr[Sf{Q)S'^{Q)p{Q 



Tr[S^{Q)S^{Q)p{Q 



Tr[^,^(0)Sf(0)p(0 



Tr[4^(0)5f(0)p(0 



Tr[4^(0)5f(0)p(0 



TT[Sf{f))S^{f))p{f) 



Pss ~l~ Pxx Pyy ~l~ Pzz')i 
Pss ~\~ pxx ~\~ Pyy Pzz')^ 
Pzs Pyx ~l~ Pxy Psz^i 
Pys ~l~ Pzx Psy ~l~ Pxz^i 
Pzs Pyx ~\~ Pxy ~\~ Psz)y 
Pxs ~l~ Psx Pzy ~l~ Pyz^ y 
Pys ~l~ Pzx ~l~ Psy ~l~ Pxz^i 
Pxs Psx Pzy ~\~ Pyz^- 



(C.5) 
(C.6) 
(C.7) 
(C.8) 
(C.9) 
(C.IO) 
(C.ll) 
(C.12) 
(C.13) 
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